Generating multi-chain configurations of an inhomogeneous melt from the knowledge 

of single-chain properties 
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Mean-field techniques provide a rather accurate description of single-chain conformations in spa- 
tially inhomogeneous polymer systems containing interfaces or surfaces. Intermolecular correlations, 
however, are not described by the mean-field approach and information about the distribution of 
distance between different molecules is lost. Based on the knowledge of the exact equilibrium single- 
chain properties in contact with solid substrates, we generate multi-chain configurations that serve 
as nearly equilibrated starting configurations for molecular dynamics simulations by utilizing the 
packing algorithm of Auhl and co-workers [J. Chem. Phys. 119, 12718 (2003)] for spatially inhomo- 
geneous systems, i.e., a thin polymer film confined between two solid substrates. The single-chain 
conformations are packed into the thin film conserving the single-chain properties and simultaneously 
minimizing local fluctuations of the density. The extent to which enforcing the near-incompressibility 
of a dense polymer liquid during the packing process is able to re-establish intermolecular correla- 
tions is investigated by monitoring intermolecular correlation functions and the structure function 
of density fluctuations as a function of the distance from the confining solid substrates. 



I. INTRODUCTION 

Many modern polymer materials consist of multiple 
components. 1 '- While the material appears homogeneous 
on a macroscopic length scale, the different constituents 
often do not mix microscopically; domains of micrometer 
size form and they are separated by internal interfaces. 
The properties of these materials arise from a complex 
interplay between the characteristics of the individual 
components in the bulk and interfacial properties. 5,4 

Therefore, abiding interest has been directed towards 
investigating inhomogeneous polymer systems, e.g., sur- 
faces of polymeric materials and polymer-solid con- 
tacts, interfaces in homopolymer blends, or microphase- 
separated morphologies in block copolymers. Computer 
simulations are an ideal tool to provide simultaneous in- 
sights into (i) the molecular conformations of the long 
flexible macromolecules at interfaces or surfaces, (ii) ther- 
modynamic or mechanical properties (e.g., interfacial 
tension) and (iii) the single-chain dynamics in the in- 
homogeneous environment and the kinetics of structure 
formation. 

The simulation of large systems of long chain 
molecules, however, poses computational challenges and 
alternative approaches, like the self-consistent field the- 
ory or polymer density functional theory, have been 
successfully applied to study inhomogeneous polymer 
systems £& 7 -' 8 -' 9 - ' 10 i 11 ' 12 i 13 ' 14 These mean-field techniques 
approximate the problem of interacting chains by the 
much simpler one of a single molecule in an external 
field that, in turn, mimics the interactions with its sur- 
rounding. Thus, significantly larger systems of longer 
chains can be described by mean-field theory than by 
computer simulations. In dense systems of long flexible 
molecules the mean-field approach yields a rather accu- 
rate description of long wave length properties?*^ The 
accuracy has been demonstrated by mapping the repre- 
sentations used in mean-field theories and computer sim- 



ulations and quantitatively comparing the results of both 
techniques ! 16 ' 17 ! 18 ' 19 ! 20 In particular, the mean-field the- 
ory provides an excellent, quantitative description of the 
molecular conformations in an inhomogeneous system. 

Given an accurate mapping between the representa- 
tions used in mean-field theories and computer simula- 
tions, it is tempting to utilize the information gained by 
mean-field theory to generate large, equilibrated multi- 
chain conformations that can be studied by computer 
simulations (e.g., to investigate the local dynamics of the 
molecules at an interface) Q However, there is a rather 
fundamental problem: The mean-field theory ignores all 
intermolecular correlations, i.e. approximating the multi- 
chain system by that of a single chain in an external 
field, all information about the relative distances between 
molecules is lost. 

In order to illustrate this problem let us consider 
the simpler case of a homogeneous dense polymer melt 
where the predictions of the mean-field theory are triv- 
ial. Within mean-field theory the chain conformations 
are unperturbed. If one generated single-chain confor- 
mations according to the pre-defined bulk distribution 
and placed them randomly into a simulation cell, the 
resulting multi-chain configuration would not represent 
a dense melt at all. The collective structure function of 
density correlations would be identical to the single-chain 
structure factor and the compressibility would be that 
of an ideal gas of chains rather than adopting the low 
value characteristic for a dense polymer liquid^ Like- 
wise, the correlations between segments belonging to dif- 
ferent molecules would be trivial instead of exhibiting a 
correlation hole that is the hallmark of intermolecular 
pair correlations in a meltf 1 ^ Since the spatial extent of 
the correlation hole is long-ranged, the chains have to 
diffuse a distance of the order of their end-to-end dis- 
tance, i? e , in order to build-up the intermolecular cor- 
relations. The necessary time, tr, increases with chain 
length and much of the benefit of using the information of 
the mean-field theory would be lost. Moreover, the corre- 
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lation hole of the intermolecular correlation function also 
gives rise to important, long-ranged intramolecular cor- 
relations in semi-dilute solutions 2 ^ and melts^ 3 - that can 
be observed in the power-law decay of the bond-bond 
correlation function along the chain contou r 24 ' 25 i 26 or a 
hump in the Kratky-plateau of the single-chain structure 
facto n 22 ' 27 It is important to emphasize that both, intra- 
and intermolecular correlations on the scale R e do not ex- 
plicitly depend on the strength of the excluded volume 
interactions or the details of the interaction potentials 
but are determined by the polymer density and chain 
cxtcnsioni 24 i 25 ' 25 



Problems due to randomly placing Gaussian chains in 
a simulation box have already been observed i 28 i 29 It was 
found that after randomly placing chains in a simula- 
tion cell and equilibrating the resulting multi-chain con- 
figuration by molecular dynamics simulations, the chain 
conformations initially swell and a time of the order of 
the Rouse time, tr, is required to achieve equilibrium. 
Furthermore, Auhl et al. demonstrated that if the chains 
are placed in the simulation box such that the density is 
uniform, much of the swelling of the conformations will 
be avoided and equilibrium will be achieved faster^ 



In the present work, we test the hypothesis that the 
incompressibility constraint, i.e. the strong reduction of 
long wave length density fluctuations in a dense poly- 
meric fluid, is sufficient to generate the intermolecular 
correlations that are ignored by mean-field theory. To 
this end, we adopt the packing algorithm of Auhl and 
co-workers 2 ^ for homogeneous systems and apply it to 
a thin polymer film confined between two structureless 
solid substrates. This simple spatially inhomogeneous 
system serves as a testing bed because it has been stud- 
ied by a variety of self-consistent field and polymer den- 
sity functional theory approachesj 19 ' 20 ^^ 3 ^' 3 ^!' 3 ^ The 
method can be readily generalized to other types of in- 
homogeneities. Moreover, rather than using Gaussian 
chain conformations, we employ single-chain conforma- 
tions that are extracted from an equilibrated melt i 19 ' 20 ' 36 
These conformations are not strictly Gaussian on large 
length scales but they already include the non-trivial 
long-ranged intramolecular correlations. 



The manuscript is arranged as follows: In the next sec- 
tion HH we describe the coarse-grained polymer model, 
the simulation techniques and the packing algorithm. In 
Sec. Mil we investigate how close the so-generated multi- 
chain configurations of a thin polymer film are to equilib- 
rium. Single-chain conformations, intra- and intermolec- 
ular correlations as well as density fluctuations will be 
considered. The paper concludes with a brief summary 
and outlook. 



II. MODEL AND TECHNIQUE 
A. Model 

We consider a widely utilized coarse-grained bead- 
spring model for polymers ; 37 ' 38 ' 39 which has successfully 
been applied to many different physical situations and for 
which many data are readily available for reference. The 
generic coarse-grained model only includes the relevant 
interactions that are necessary to describe the universal 
behavior of long, flexible macromolecules in a dense solu- 
tion. Harsh repulsive excluded volume interactions and 
longer-ranged, attractive van-der-Waals attractions act 
between all segments of all chains and are modeled by 
a spherically truncated and shifted Lennard- Jones (LJ) 
potential, 



U hJ (r) 
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e and a define the units of energy and length, respec- 
tively. The cutoff radius, r c = 2 x 2 1 / 6 <r, is given by 
twice the position of the minimum in Eq. (H]). Addition- 
ally, the connectivity along the macromolecule is modeled 
by a finitely extensible nonlinear elastic (FENE) poten- 
tial that acts between bonded neighbor segments along 
the same polymer chain, 
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where i?o = 1.5 a is the maximum extension of a bond 
and k = 30 e/ a 2 is the spring constant. 

Polymer segments interact with the unstructured solid 
substrate via an integrated Lennard- Jones 3-9-potential 



U wa u(z) — ££w 



(3) 



where the dimensionless Hamaker constant £w = 3.0 is 
a measure for the attraction between solid substrate and 
polymer melt. 

Immediately after rearranging the chains by means of 
the packing algorithm, which will be described below, 
there still remains some overlap between different seg- 
ments. In this case, the use of Eq. (pj for non-bonded 
interactions leads to numerical instabilities. At this ini- 
tial stage, we therefore employ for all non-bonded interac- 
tions a potential with a maximum force at short distances 
r < r ic . 



\U-Lj[r), r>r fc 

In the limit of vanishing rf c — > 0, the potential converges 
to the full interaction potential^ 



3 





N 


n 


L z /a 




Dt/<t 2 


Ri/a 2 


R g /a 


tr/t 


I 


64 


300 


31.56 


0.611 


4.36 • 10 _a 


104 


4.19 


805 


II 


128 


230 


34.62 


0.710 


1.27 • 10~ 3 


209 


5.88 


5558 


III 


256 


230 


43.61 


0.710 


2(1) • 10" 4 


446 


8.74 


8(4)-10 4 



TABLE I: Parameters of the simulations. N: number of ef- 
fective segments (interaction centers) per chain. R e : mean- 
squared end-to-end distance. D: self-diffusion coefficient of 
the molecules' center of mass, tr: Rouse times calculated via 
t r = Ri/(3n 2 D)^ 



In the following, we study a monodisperse melt of n 
polymers that are comprised of N effective segments per 
chain in the canonical ensemble. We use a temperature of 
ksT/e = 1.68 where fcg denotes Boltzmann's constant. 
The melt is confined between two impenetrable unstruc- 
tured solid substrates parallel to the xy-plane. The solid 
substrates are located at z = —L z /2 and z — L z /2. Pe- 
riodic boundary conditions are applied in the two other 
directions, x and y. The mean segment number density, 
p, as well as the other parameters are compiled in table 

in 

Two different methods have been utilized to equilibrate 
the configurations: On the one hand, a Monte Carlo 
(MC) program with a mixture of both, slithering snake 
and random local moves, is employed because these MC 
moves can quickly equilibrate even long chains at moder- 
ate density. On the other hand, we use a parallel Molecu- 
lar Dynamics (MD) program with Nose-Hoover thermo- 
stat to simulate at constant temperature. 41 ! 42 Within this 
technique, the system is coupled to a heat bath, which is 
represented by one additional degree of freedom, 77, with 
canonically conjugated momentum, p^. The concomitant 
equations of motion are 

i-i = —, p i = Fi - ^pi, (5a) 
m Q 

nN 2 

? ?=77' Pn=J2 — -9 k BT. (5b) 
Q ^ m 

i] fluctuates around zero; its "mass", Q, is chosen as to 
realize a canonical ensemble, to is the mass of a segment 
and we set to = 1 in the remainder of this article, g — 
3nN denotes the number of degrees of freedom. The 
equations of motion {5| were integrated using the velocity 
Verlet algorith m 43 ! 44 with a time step of At = 0.002 r, 
where r = aymje denotes the time unit in Lennard- 
Jones parameters. 



B. Packing process 

In this section, we describe our packing algorithm, 
which we have used to create polymer melts of long chains 
that are already very close to equilibrium. This algorithm 
is based on an idea of Auhl and co-workers^ Here, we 
generalize it to spatially inhomogeneous systems. The 



underlying hypothesis is to create intermolecular corre- 
lations by enforcing incompressibility of the dense poly- 
mer solution. The strategy of the packing algorithm con- 
sists in reducing local fluctuations of the segment den- 
sity without corrupting the intramolecular correlations. 
This algorithm has been implemented as a Monte Carlo 
program performing trial moves on a priori equilibrated 
chain conformations, which are treated as rigid bodies. 

The equilibrated polymer melt in contact with the solid 
substrate exhibits ^-dependent intra- and intermolecular 
correlations on different length scales i 45 ' 46 The spatial in- 
homogeneity makes correlation effects also visible in the 
spatial dependence of single-chain properties (e.g., the 
segment density distribution). On the one hand, there 
are intricate packing effects of the fluid of segments that 
manifest themselves in oscillations of the density profile 
as a function of the distance, z, from the solid substrate. 
In a dense melt, this local arrangement of segments de- 
cays on the length scale of the segment size or the screen- 
ing length, £, of density fluctuations, which is micro- 
scopic. Our packing algorithm does not aim at generating 
this short-ranged fluid ordering because this local struc- 
ture on the segmental scale equilibrates on time scales 
that are fast compared to the Rouse time of the macro- 
molecules. On the other hand, the presence of the solid 
substrate influences the polymer structure over distances 
up to the molecular size. In this region - denoted inter- 
phase - the polymer chains align with the solid substrate, 
the back-folding of the conformations deepens the corre- 
lation hole of the intermolecular pair correlation func- 
tion and the density of the molecules' center-of-masses 
depends non-trivially on the distance from the solid sub- 
strate. These effects have to be captured by the pack- 
ing algorithm because their equilibration requires the 
molecules to diffuse a distance of the order of the molec- 
ular extension and the concomitant time scale is set by 
the Rouse-time.— 

Our packing algorithm requires (i) equilibrated single- 
chain conformations of a bulk system that do incorporate 
the non-trivial long-ranged intramolecular correlations 
imparted by the incompressibility constraint and (ii) the 
distribution of the single-chain conformations in the spa- 
tially inhomogeneous system. The latter property can 
be obtained with high accuracy by self-consistent field 
or polymer density functional calculations evaluating the 
single-chain partition function via a partial enumeration 
over a large ensemble of bulk single-chain conformations. 

In order to assess the quality of the packing algorithm, 
however, we adopt in the following a different procedure 
that does not invoke the additional approximations of the 
self-consistent field or polymer density functional calcu- 
lations. We equilibrate a dense polymer melt in the pres- 
ence of the solid substrate. This equilibrated multi-chain 
configuration serves (i) as a reference to which we com- 
pare the result of the packing process and (ii) it provides 
exact information about the distribution of single-chain 
conformations in the spatially inhomogeneous system. 
We extract single-chain conformations at a given distance 
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from the spatial inhomogeneity from equilibrated melt 
configurations (at different times) and our packing algo- 
rithm re-assembles them to a multi-chain configuration. 
Therefore, the deviations of the multi-chain configuration 
from the equilibrated melt solely stem from the packing 
algorithm. 

In the first stage, we generated an initial multi-chain 
configuration by randomly drawing equilibrated molecu- 
lar conformations from a large ensemble of single-chain 
conformations. The distance from the spatial inhomo- 
geneity (i.e., the z-coordinate) and the orientation of each 
single-chain conformation in the newly assembled config- 
uration is identical to the value in the equilibrated con- 
figuration and the density profiles of both configurations 
are also statistically identical. We emphasize that only 
information about single-chain properties, which can be 
obtained from a mean-field theory enter the packing pro- 
cess. This first stage generates an initial configuration 
that reproduces the single-chain properties of the equili- 
brated melt (e.g., the molecules' center-of-mass distribu- 
tion) but the intermolecular correlations parallel to the 
solid substrate are that of an ideal gas of polymers. 

In the following packing procedure, single-chain confor- 
mations are treated as rigid bodies (i.e., the intramolec- 
ular correlations are frozen) and the distance to the 
(nearest) solid substrate as well as the orientation of the 
molecule is kept fixed. In order to introduce intermolec- 
ular correlations in the direction parallel to the spatial 
inhomogeneity, we minimize local density fluctuations. 
To this end, we define a local number of neighbors, rii >s 
of the s th segment of polymer i by counting the number 
of segments within a sphere of radius r p around the refer- 
ence segment^ This procedure guarantees translational 
invariance and isotropy in the bulk. In a spatially inho- 
mogeneous system, however, the local number of neigh- 
bors and its fluctuations depend on the distance, z, from 
the solid substrate. In practice, we discretize the system 
into 2nf, + 1 bins such that there are rib bins of constant 
width Afc on each side. For the multi-chain configura- 
tion, we define the local number of neighbors in layer Zk 
with &=!,••■ , 2nb + I according to 



n(zfc) 



I dz £Li Ss=i "v^fc - * 
I dz E"=i EfLi - z i) 



(6) 



where the integral runs over the fc th bin. The variance 
of the local number of neighbors is a measure of density 
fluctuations and we minimize it for each layer parallel 
to the solid substrate. Thus, the target function for the 
second step of the packing process is given by 

2n b + l 

E=Y, A kE(z k ) with E(z k ) = (n 2 )(z)-(n) 2 (z). 
k=l 

(7) 

To preserve the structure of the inhomogeneous sys- 
tem, the distance between every segment and the inho- 
mogeneity must not change during the second stage of 



the packing process. The inhomogeneity under consider- 
ation is both invariant under translations in the parallel 
directions, x and y, and rotations around the perpendic- 
ular z-axis. Additionally, it is invariant under reflections 
at the middle plane of the thin film. In order to reduce lo- 
cal density fluctuations and simultaneously preserve the 
structure in z-direction, we utilize the following set of 
moves in the second stage: (i) Translation of a single chain 
parallel to the solid substrates, (ii) Rotation of a single 
chain by an arbitrary angle along an axis perpendicular 
to the solid substrates, (iii) Reflection of a chain by the 
middle plane of the system, (iv) Replacement of a chain 
in the system by a different chain at a different position 
according to the mean-field prediction of the single-chain 
properties. Such a change of a single-chain conformation 
is accepted if it will decrease the target function, E. 

The typical values we have employed are = 5 re- 
sulting in A = 3.0 a. The second stage of the packing 
process started with a value of r p = 4.0 a, which was 
successively decreased in steps of Ar p = 1.0 a down to 
r p = 2.0(7. The typical runtime of the packing process 
for the system with n = 300 and N = 64 with a total 
number of 4 • 10 6 MC steps on a single Intel Xeon 5160 
CPU is roughly 3 h. 

After this second stage of the packing process, the 
multi-chain configuration is utilized as starting config- 
uration of a MD run. Since, immediately after the sec- 
ond stage, there still remain some local overlaps of the 
segments, we cannot use the full interaction potential 
but rather employ the linearized potential (UJ) for non- 
bonded interactions. The initial value of rf c = 0.9 a is 
gradually decreased by 0.001 a every 20 r until a final 
value of rf c = 0.82 r was reached. This "push-off" is 
very effective in removing local segmental overlaps. This 
third stage completes the generation of an overlap-free, 
inhomogeneous multi-chain configuration for our coarse- 
grained bead-spring model. 

In order to assess how close to equilibrium the so- 
generated configuration is, we studied the configurations 
with the full interactions between all segments for ad- 
ditional 3400 t by MD simulations. Initially, we as- 
signed random velocities to each segment according to 
the Maxwell-Boltzmann statistics such that the temper- 
ature in the initial state already matched the target tem- 
perature. 



III. RESULTS AND DISCUSSION 



A. Packing Process 



We apply the packing process described in Sec. Ill Bl 
to confined polymer melts, which are comprised of 
molecules with different chain lengths, N — 64, 128, 256. 
In this subsection, we describe the impact of this pack- 
ing scheme on the properties of the polymer melt with 
N = 64 by inspecting the correlation hole and the collec- 
tive structure function. 
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The collective structure function of the density is given 



G{q\\) = 



nN 



n N 

EE< 

»=1 s=l 



(8) 



It is a measure of the strength of density fluctuations on 
all length scales. In Eq. ((5]), qn denotes a wave vector par- 
allel to the solid substrate and q\\ = |q|| |. Due to the pe- 
riodic boundary conditions in the two parallel directions, 
the components of are restricted to q^j = 2irn i /L i 
with m = 0, 1, 2, ... . Wave vectors are measured in units 
of the molecules' radius of gyration, R g . 

In Fig. [TJ we present G(q\\ ) during the second stage 
of the packing process with r p — 4.0 ct, 3. Oct, and 2.0 ct, 
respectively. The data after the first stage of the packing 
process and in equilibrium are also displayed for com- 
parison. The values of q\\ corresponding to 2r p (i.e., the 
diameters of the sphere around a segment) are indicated 
by vertical dashed lines in the figure. 

Since the lateral arrangement of single-chain confor- 
mations is random after the first stage of the packing 
process, G(q\\ ) is exactly given by the structure function 
of a single chain 




N 

E< 

s=l 



(9) 



The peak at large wave vectors, q\\R g ~ 32, is a con- 
sequence of the intramolecular packing and reflects the 
length scale of a single bond. In an equilibrated dense 
melt, this peak also exists but it stems from both, intra- 
and intermolecular packing. In our model, there is a 
small mismatch between the preferred bonded and non- 
bonded distances^ which explains the small difference 
between this local packing peak after the first stage of 
the packing and in an equilibrated melt. 

During the second stage of the packing process, den- 
sity fluctuations on the length scale r p are reduced and 
r p is successively decreased. The impact of the packing 
algorithm can be observed quite clearly by the presence 
and shift of the additional peak in G(q\\) at intermedi- 
ate wave vectors. The packing reduces density fluctua- 
tions for <C 2n/r p but leaves fluctuations on shorter 
length scales as well as the single-chain structure unal- 
tered. Thus, the collective structure function develops 
a peak around gjj rs 27r/2r p , which shifts to larger wave 
vectors as we decrease r p . 

After packing with the largest value, r p = 4.0 ct, the 
density fluctuations have been significantly reduced on 
length scales q\\R g < 4; in fact, they are even smaller 
than in equilibrium. In the course of the subsequent 
packing, the fluctuations for small gii slightly increase 
again because both long-ranged and shorter-ranged fluc- 
tuations are simultaneously minimized. After packing 
with r p = 3.0ct the peak in G has moved to q\\R g ~ 5.8 of R g < z < 2R g and the last one in the center is rep- 
and its intensity has further decreased. After the final resentative for the bulk behavior. In Fig. [2] we show the 



FIG. 1: Comparison of the collective structure function, G(q«) 
(TV = 64, n = 300, pa i = 0.611), of the system in equilibrium 
with the packed system before, after the completion of each 
packing step with a certain radius r p , and after application of 
mixed moves. S(q) denotes the single-chain structure factor. 



packing with r p = 2.0 ct, this peak decreases still further, 
however, the small length scales remain unaffected by 
the packing process, i.e. G(q\\) still resembles the single- 
chain structure factor, S(q), for q\\R g > 8. Any further 
decrease of r p does not improve the situation. Instead of 
moving G(qu) closer to its equilibrium value, we observe 
that long wave length density fluctuations increase, i.e., 
the small q-regime is getting worse in favor of the large 
g-regime. Since the relaxation times for large q are small, 
we decide to terminate the second stage of the packing 
process at r p = 2.0 ct. 

A comparison of these structure functions with G(q\\ ) 
of an equilibrated melt (cf. Fig. |lj) demonstrates that the 
second stage of the packing process is able to decrease 
the long wave length fluctuations significantly and that 
the collective structure function for small values of qu 
approaches the equilibrium value. 

In passing, we note that the limiting value G(q\\ — *• 0) 
is not fixed by the packing process, but can be adjusted 
to values between that of the second stage of the pacing 
process at r p = 2.0 ct and that of r p = 4.0 ct by mixing 
moves with different values of r p . Choosing a ratio of 9 : 1 
of moves with r p = 2.0 ct and 4.0 ct, we observe in Fig. [T] 
that the position of the packing peak levels off, however 
the peak broadens slightly towards smaller values of q\\ . 

We also investigate the packing as a function of the dis- 
tance from the spatial inhomogeneity. Since the relevant 
length scale for local packing is given by the correlation 
length £ (£ » 2.5 a for N = 64) and the length scale for 
intermolecular packing is given by the radius of gyration, 
R g , we divided the system into four slabs: The first in- 
terval from < z < £ is close to the solid substrates, the 
second spans £ < z < R g , the third region is comprised 
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FIG. 2: G(gy) after packing with r p — 2.0 a for the intervals 
specified in the text (N = 64). The noticeable discrepancy 
between the equilibrium distribution and that from the pack- 
ing algorithm at < z < £, poses no severe problems, since 
the distribution of the molecules center-of-mass is correct by 
construction and the segments will move quickly into this vol- 
ume. 



collective structure function, G(q\\), in these four inter- 
vals after the second stage of the packing process. We 
observe that G(q\\ ) in the interval from < z < £ devi- 
ates from its equilibrium distribution. This deviation is 
expected because our algorithm cannot reproduce local 
packing on the segmental scale. The closer the slab is to 
the solid substrate, the larger are fluctuations at small 
qu. For distances z > 2R g , the fluctuations resemble the 
structure function of Fig. [T] that has been obtained by 
averaging over the entire sample thickness. 

While the structure function measures the density fluc- 
tuations on different length scales it does not distinguish 
between intermolecular packing and intramolecular con- 
formations. To quantify the distribution of distances be- 
tween different chains, we calculated the intermolecular 
radial pair correlation function for segments within the 
same distance Az from the solid substrate 



n N 



Pointer ( 



i^j s,t—l 



. (10) 

as well as the total radial pair correlation function of 
the segments within the same distance from the solid 
substrate 

^ n N 

p9 ^0 = 7^f 51 ( s - r Jt\ - r\i) 6 {z is - z jt )) . 



nN 

s,t— 1 

(ii) 

The former exhibits the intermolecular correlation hole 
while the latter is the Fourier transform of the struc- 
ture function, G. In the main panel of Fig. [3J we 
present Winter (f\\) averaged over the entire system. Af- 
ter the first stage of the packing process (initial config- 



uration), there are no intermolecular correlations, i.e., 
Sinter M = 1- Already after the second stage of the 
packing with r p = 4.0 a, a pronounced correlation hole 
is established, and for rii > 5c, Sinter(?~||) matches the 
equilibrium distribution already very well. Upon de- 
creasing the value of r p the correlation hole is getting 
deeper and Sinter (^ii) is approaching its equilibrium dis- 
tribution. Two different regions can be distinguished: (i) 
F° r r \\ ^ 5 a, the intermolecular correlations are deter- 
mined by the packing of whole molecules and this effect is 
correctly reproduced by the packing algorithm. Hence, 
deviations from the equilibrium distribution are small. 
This is an important regime because creating this corre- 
lation hole on the scale of R g by a diffusive motion of the 
molecules would require an equilibration time of the or- 
der tr. (ii) For rii <C 5cr, the intermolecular correlations 
exhibit local packing which gives rise to an oscillatory 
behavior. For r« < 1 a, these local packing effects domi- 
nate the behavior and the Sinter (rii ) after the second stage 
deviates significantly from the equilibrium distribution. 
The length scale of the local packing effects is set by 
the size of a segment and the corresponding relaxation 
times are fast. For the rather small number of segments, 
N = 64, the two regimes are not clearly separated and 
there is a rather gradual cross-over between the limiting 
behaviors that characterize the local structure of the fluid 
of segments and the polymeric correlation hole. 

The inset of Fig. [3] depicts the total pair correlation 
function, <7(rii). At the beginning of the second stage 
of the packing process, g(r\i) exceeds unity for distances 
of the order R g , i.e., <?(rii) — 1 is proportional to the 
intramolecular correlations. After the packing, however, 
g(r\\) has decayed to unity beyond r\\ — ha w 2£ and 
the whole distribution closely resembles the equilibrium 
distribution. 

In Fig.[4]we analyze the spatial correlations, <7i n ter(f||), 
in slabs with different distances from the solid substrate. 
The slabs are defined in the same way as in the discussion 
of the collective structure function. Although Sinter 
shows in all regions a qualitatively similar behavior, the 
correlation hole near the solid substrates is considerably 
deeper and longer-ranged than in the bulk because of 
the back-folding of the chain conformations. On large 
scales (rii > 4er), this effect is perfectly reproduced by 
our packing algorithm. In the inset of Fig. 31 we present 
7~]j • [l — Sinter (?"|| )] to highlight the dependence of the 
correlation hole on the distance to the solid substrate. 



B. Fast relaxation 

After having generated multi-chain configuration with 
the first two stages of our packing algorithm, we per- 
formed MD simulations to investigate the final equilibra- 
tion during which the local packing structure of the fluid 
of segments is established. Due to the remaining over- 
lap of segments, we use a push-off procedure (cf. Sec. |TTJ) 
with the linearized Lennard- Jones potential of Eq. (U) for 
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FIG. 3: The intermolecular pair correlation function (/inter (f y) 
after each packing step for the system = 64 in comparison 
with the equilibrium distribution. The depth of the inter- 
molecular part of the correlation hole is increasing with every 
packing step and on lengths larger than the rii > 4cr the 
final distribution matches the equilibrium distribution very 
well. In the inset we present the total pair correlation func- 
tion <7(V||) including the intramolecular contributions. 




FIG. 4: (/inter for N = 64 after packing with r v — 2.0 a 
for intervals with a different distance to the solid substrate. 
Sinter (r || ) shows in the interval < z < £ a deeper correlation 
hole, that is reproduced by the packing algorithm. The inset 
shows r|j • [l — <?intcr(V||)] so that the correlation hole becomes 
clearly visible. 



non-bonded interactions. During these simulations, we 
monitored the collective structure function, G(q\\), the 
intermolecular pair correlation function, <?inter(?"||)7 and 
the mean-squared size (R 2 (s)) of chain segments of curvi- 
linear length s = \n — m\ < N. Note that R 2 includes 
the parallel and perpendicular extension of the molecules 
in confinement. The latter quantity is averaged over all 
segments of size s, where n and m denote segment in- 
dices along the chain. Since (R 2 (s)) is given by the in- 
tegrated bond-bond-correlation function, it serves as a 




1000 



FIG. 5: {R 2 (s))/s = (R 2 {s) + R 2 ± (s))/s, mean-squared size 
of chain segments of curvilinear length s normalized to s in 
equilibrium. 



sensitive measure for intramolecular correlations.— For a 
Gaussian chain (R 2 (s))/s = b 2 , where b is the statistical 
segment length. For chains in a melt, there are devia- 
tions from this law, particularly at small s. In Fig. we 
have plotted (R 2 (s))/s in equilibrium for the three chain 
lengths under consideration. For N = 64 and N = 128, 
both distributions show similar behavior and they reach 
a plateau at b 2 for s — > N. The distribution for N = 256 
shows the plateau at a higher value of b 2 , i.e. the chains 
are slightly more swollen because of the confinement. 

In our MD runs we observed that the relaxation on 
length scales of the size of a few bond lengths occurs 
very rapidly. Inspecting G(q\\), we notice that already 
after At = 10 r there is hardly any difference compared 
to the equilibrium distribution visible. Even density fluc- 
tuations around q\\R g ~ 7 equilibrated completely, and 
also g(ru) shows the same fast relaxation to equilibrium. 
In Fig. [B] we present the fast equilibration of G(q\\) for 
the slab < z < £, which has reached approximately the 
equilibrium distribution after 0.1 tr. 

Regarding (R 2 (s)), however, it turns out that the in- 
tramolecular correlations have been disturbed and that 
their relaxation takes considerably longer than the re- 
laxation of the local packing. The main panel of Fig. 
shows (R 2 (s, t)) as a function oit normalized by the equi- 
librium distribution (R 2 q (s)} for N = 64. The figure 
presents data for several values of s. tr is the Rouse 
time of a chain^S (cf. Tab. Q|. The strongest devia- 
tions from equilibrium occur at s rs 15 and they are 
smaller than 8 %. Larger scales (larger values of s) show 
smaller deviations which decay exponentially in time. Af- 
ter roughly 1.5 tr, the deviations are smaller than 1 % on 
all length scales and we consider the intramolecular cor- 
relations as being equilibrated. The inset of Fig.[7]depicts 
(R 2 (s , t)) I (R 2 q {s)) as a function of s for fixed values of 
t. The peak of the strongest deviations moves slowly to 
larger s, importantly, it does however not propagate to 
the largest length scales. 
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FIG. 6: G(qu) for the slab < z < £ at the beginning of the 
MD runs (N = 64). 



1.08 



1.06 



04 

V 



1.04 



1.02- 




FIG. 7: The main panel presents (R 2 (s, t))/{R 2 q (s)) for N = 
64 as a function of time for a few fixed values of s, while in 
the inset the same quantity is depicted as a function of s for 
various time steps. 



In Fig. [H we finally compare the temporal evolution 
of (R 2 (s,t))/{R^ q (s)) after the first stage of the packing 
with a configuration generated after the second stage. 
With a maximal deviation from the equilibrium distri- 
bution of 26% at s ea 13 after the first stage, the in- 
tramolecular correlations are much stronger perturbed 
and also occur on slightly larger values of s than in the 
configuration after the second stage. But even worse, 
these deviations in the configuration after the first stage 
propagate to larger length scales and therefore increase 
the required relaxation time significantly. 



IV. CONCLUSIONS AND OUTLOOK 

We have investigated the properties of a packing algo- 
rithm based on an idea by Auhl and co-workers^ that 
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FIG. 8: Comparison of the temporal evolution of 
(R 2 (s,t))/(Rc q (s)) between an initial configuration gener- 
ated by the packing algorithm (left) and a randomly packed 
configuration (first stage, right) for some fixed time steps t 
(N = 128). 



arranges given single-chain conformations extracted from 
an inhomogeneous polymer melt in order to generate a 
multi-chain configuration, which can serve as a start- 
ing configuration for a molecular dynamics simulation. 
The data demonstrate that reducing local density fluc- 
tuations, which are quantified by the average number of 
segments that surround a reference bead, we are able 
to re-establish intermolecular correlations on the length 
scale of the molecules' extension, R e . Most prominently, 
the packing algorithm generates the correlation hole in 
the intermolecular pair correlation function. 

In order to assess the quality of the packing algo- 
rithm, we utilize single-chain conformations extracted 
from an equilibrated thin polymer film. The so-generated 
single-chain properties are exact by construction. In 
practice, mean-field techniques that evaluate the single- 
chain partition function via a partial enumeration over a 
large ensemble of equilibrated single-chain conformations 
of the bulk will provide an excellent approximation for 
the single-chain properties in an inhomogeneous system. 
Since the single-chain conformations are extracted from 
the bulk they already incorporate the long-ranged in- 
tramolecular correlations that the incompressibility con- 
straint imparts on the conformations (e.g., the power-law 
decay of the bond-bond correlation function along the 
molecule). 

The local fluid-like packing structure on the length 
scale of the excluded volume screening length is not re- 
produced by the packing algorithm. This local structure, 
however, can be established on a short time scale com- 
pared to the relaxation time of the inter- and intramolec- 
ular correlations on the scale R e . Thus, the algorithm 
is particularly useful for creating initial configurations 
for particle-based simulations of inhomogeneous polymer 
melts or semi-dilute solutions on the basis of mean-field 
results. 
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There remain two problems to be addressed: (i) For 
inhomogeneous melts of very long chains, even the gen- 
eration of equilibrated bulk configurations is a com- 
putational challenged Gaussian chain conformations, 
however, are inappropriate because the long-ranged in- 
tramolecular correlations require protracted long times 
to equilibrate. Instead of drawing long polymers di- 
rectly from an equilibrated melt it might be feasible 
to generate appropriate conformations by joining short 
chain o 47 ' 48 accounting for the weak effective intramolec- 
ular interactions 4 ^ due to the incompressibility. This pro- 
posed procedure is very similar to the analytical pertur- 
bation calculations performed by Semenov and Johncr 23 
and we expect it to reproduce the non-trivial long-ranged 
intramolecular correlations.— These single-chain confor- 
mations can be used as a starting point for mean-field 
calculations of inhomogeneous systems and the results of 
the latter serve as input into the packing algorithm. 

(ii) We have applied the packing algorithm to a poly- 
mer melt with the simplest inhomogeneity, i.e. a thin 
polymer film confined by symmetric solid substrates. We 



emphasize that the packing procedure can be readily gen- 
eralized to more complicated geometries. The moves, 
which are utilized to arrange the single-chain confor- 
mations during the packing process, must not alter the 
single-chain properties (i.e., the intramolecular structure 
and the position and orientation with respect to the spa- 
tial inhomogeneity). Thus, in the case, where transla- 
tional invariance is broken in all three spatial directions, 
only replacement moves are acceptable but no transla- 
tions or rotations. 
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